1, 2012 0:10 WSPC/INSTRUCTION FILE nishimura7 



International Journal of Modern Physics: Conference Series 
© World Scientific Publishing Company 



(3+l)-DIMENSIONAL EXPANDING UNIVERSE FROM A 
LORENTZIAN MATRIX MODEL FOR SUPERSTRING THEORY 

IN (9+l)-DIMENSIONS 



JUN NISHIMURA 

KEK Theory Center, High Energy Accelerator Research Organization 
1-1 Oho, Tsukuba, Ibaraki 305-0801, Japan 
jnishi@post.kek.jp 

We study the Lorentzian version of the type IIB matrix model as a nonperturbative 
formulation of superstring theory in (9+l)-dimensions. Monte Carlo results show that 
not only space but also time emerges dynamically in this model. Furthermore, the real- 
time dynamics extracted from the matrices turns out to be remarkable: 3 out of 9 spatial 
directions start to expand at some critical time. This can be interpreted as the birth of 
our Universe. 
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1. Introduction 

Fundamental questions concerning our Universe are: 

(i) Why does our Universe have (3+l)-dimensions? 

(ii) Why is it expanding ? 

We provide explicit answers to these questions from a nonperturbative formulation 
of superstring theory in (9+l)-dimensions pQ. 

More specifically, we use the type IIB matrix model [2J, which was proposed as 
a nonperturbative formulation of string theory based on type IIB superstrings in 
10 dimensions. The Euclidean version of this model has been studied by various 
approaches. In particular, the Gaussian expansion method proposed in Ref. [3] was 
used to calculate the free energy of the SO(d) symmetric vacua (2 < d < 7) and it 
was found that d — 3 gives the minimum [4]. This implies that the SO(10) symmetry 
of the Euclidean model is spontaneously broken down to SO(3). Moreover, the extent 
of space-time in the extended d directions and that in the shrunken (10— d) directions 
turn out to have a finite ratio even in the large- N limit [4] . While these results reveal 
interesting dynamical properties of the Euclidean model, the connection to our real 
space-time is not very clear. This motivated us to consider the Lorentzian model pQ. 

Let us quote here an interesting statement made by Seiberg in his rapporteur 
talk entitled "Emergent Spacetime" [5] at the 23rd Solvay Conference in Physics 
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in 2005: Understanding how time emerges will undoubtedly shed new light on some 
of the most important questions in theoretical physics including the origin of the 
Universe. Indeed, we will see that the emergent spacetime is naturally realized in 
the Lorentzian matrix model, and that the results can be interpreted as describing 
the birth of our Universe. 

The rest of this article is organized as follows. In section [2] we review previous 
works on type IIB matrix model. In section[3]we define the Lorentzian matrix model. 
In section |4] we present our Monte Carlo results for the Lorentzian matrix model. 
Section [5] is devoted to a summary and discussions. 

2. Previous works on type IIB matrix model 

The action of the type IIB matrix model is given by [2] 



where (fj, = 0, • • • ,9) and ^> a (a = 1, • • • , 16) are NxN Hermitian matrices. The 
Lorentz indices /i and v are contracted using the metric rj — diag(— 1, 1, • • • ,1). The 
16 x 16 matrices L M are ten-dimensional gamma matrices after the Weyl projection, 
and the unitary matrix C is the charge conjugation matrix. The action has manifest 
SO(9,l) symmetry, where and ty a transform as a vector and a Majorana-Weyl 
spinor, respectively. 

There are various evidences that the model gives a nonperturbative formulation 
of superstring theory First of all, the action ([1} can be viewed as a matrix 
regularization of the worldsheet action 



of type IIB superstring theory in a particular gauge known as the Schild gauge. 
(The "Poisson brackets" are defined here by {/(£)>.9(£)} = e ij§£"§§~-) ^ has also 
been argued that configurations of block-diagonal matrices correspond to a collec- 
tion of disconnected worldsheets with arbitrary genus. Therefore, instead of being 
equivalent just to the worldsheet theory, the large- N limit of the matrix model is 
expected to be a second-quantized theory of type IIB superstrings, which includes 
multi-string states. Secondly, D-branes are represented as classical solutions in the 
matrix model, and the interaction between them calculated at one loop reproduced 
correctly the known results from type IIB superstring theory |2 j . Thirdly, one can de- 
rive the light-cone string field theory for the type IIB case from the matrix model [8] 




(1) 
(2) 




Sschild 




a There are also other proposals for nonperturbative formulations of superstring/M theory based 
on matrix models such as Matrix Theory [B] and Matrix String Theory [7]. 
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with a few assumptions. In the matrix model, one can define the Wilson loops, which 
can be naturally identified with the creation and annihilation operators of strings. 
Then, from the Schwinger-Dyson equations for the Wilson loops, one can actually 
obtain the string field Hamiltonian. 

In all these connections to string theory, it is crucial that the model has two 
kinds of fermionic symmetries given by 

5^A a = ^ el C^ A1 * , 

<5 (1) * =iT^[A ll ,A v ]e l , W 
S^A, = , 

5^ =e 2 l, W 
where 11 is the unit matrix. It also has the bosonic symmetry given by 

(6^A, = c,l, 

{s^v =0. [V 

Let us denote the generators of (0, ([6]) and ([7]) by Q^ 1 ', and P M , respectively, 
and define their linear combinations 

Q (1) =g (1) + g (2) , Q^=*(Q«-Q< 2 >). (8) 

Then, we find that the generators satisfy the algebra 

[eiCQ« e 2 CQ^] = -2^' ei CF^ 2 P M , (9) 

where i, j — 1,2. This is nothing but the ten-dimensional N — 2 supersymmetry. It 
is known that field theories with this symmetry necessarily include gravity, which 
suggests that so does the type IIB matrix model. When we identify © with the 
ten-dimensional M = 2 supersymmetry, the symmetry ((7|) is identified with the 
translational symmetry in ten dimensions, which implies that the eigenvalues of 
A^ should be identified with the coordinates of ten-dimensional space-time. This 
identification is consistent with the one adopted in stating the evidences listed in 
the previous paragraph, and shall be used throughout this article as well. 

An interesting feature of the type IIB matrix model is that the space-time itself 
is treated as a part of dynamical degrees of freedom in the matrices. Therefore, it is 
possible that a four-dimensional space-time is generated dynamically in this model. 
This issue has been studied in the Euclidean version of the model j3H4|[9Tl23) , which 
can be obtained from ([l} by making a Wick rotation 

A Q h> iA w , F° i ^ -iT 1Q . (10) 

Note that the Euclidean model has SO(10) symmetry instead of SO(9,l). 

In order to discuss the spontaneous symmetry breaking (SSB) of SO(10) in the 
large- N limit, we consider the "moment of inertia" tensor [91IT0] 

1 

N' 



T MV = -tr(A ll A v ), (11) 



October 1, 2012 0:10 WSPC/INSTRUCTION FILE nishimura7 



4 Jun Nishimura 




Fig. 1. (Left) The free energy of the SO(d) symmetric vacuum is plotted against d. The horizontal 
line represents the value / = log 8 — ? = 1.32944... obtained from the conjecture by Krauth, 
Nicolai and Staudacher | 24 | . (See footnote [b]) The dotted line connecting the data points is drawn 
to guide the eye. (Right) The extent of space-time B? and r 2 in the extended and shrunken 
directions, respectively, are plotted against d. The solid and dashed lines connecting the data 
points are drawn to guide the eye. 



which is a 10 x 10 real symmetric tensor. We denote its eigenvalues as Xj (j — 
1, ■ ■ • , 10) with the specific order 

Ai > A 2 > ■•• > Aio . (12) 

If the SO(10) symmetry is not spontaneously broken, the expectation values (Xj) 
(j = 1, • • • , 10) should be all equal in the large- N limit. The dynamical generation 
of (i-dimensional space-time corresponds to an SO(d) symmetric vacuum, in which 
the expectation values (Xj) (j — 1, ■ ■ • , d) are equal in the large- N limit, but the 
remaining ones (Xj) (j = d + 1, • • • , 10) are much smaller. 

Let us show recent results obtained by the Gaussian expansion method [4]. In 
Fig. Q] (Left) we plot the free energy of the SO(d) symmetric vacuum for 2 < d < 7 
at order 3 of the expansion. The result decreases monotonically as d decreases from 
7 to 3, and it becomes much larger for d — 2. Thus, the SO(3) symmetric vacuum 
gives the smallest free energy, which suggests the SSB of SO(10) down to SO(3)@ 

Let us discuss the results for the extent of space-time represented by the eigen- 
values (fl"2|) . In Fig. Q] (Right) we plot the result for the extended directions (R 2 ) 
and the shrunken directions (r 2 ) for each d. We find that r 2 stays almost constant 
at r 2 = 0.1 ~ 0.15, which seems to be universal for all the SO(d) symmetric vacua 
with 2 < d < 7. On the other hand, the results for R 2 are found to be larger for 
smaller d. Thus the extent of space-time seems to be finite in all directions. While 
the observed SSB of SO(10) is an interesting dynamical property of the Euclidean 



b The d-dependence of the free energy is quite analogous to the one observed in the six-dimensional 
model [25] . There the value of the free energy tends to decrease slightly as one goes from order 
3 to order 5. Considering such artifacts due to truncation, we speculate that the Krauth-Nicolai- 
Staudacher conjecture | 24 | actually refers to the partition function for the SO(10) symmetric 
vacuum. See Fig. [T] (Left). 
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model, the connection to the real space-time is unclear. 
3. Defining the Lorentzian matrix model 

The most crucial difference between the Euclidean and Lorentzian models is the 
bosonic part of the action ([2]l, In the Euclidean model, it becomes positive semi- 
definite 

5 b cx tr {F^f , (13) 

where we have defined Hermitian matrices F^ = i[^4 M , A u ], There is a flat direction 
corresponding to [A^, A v \ ~ 0, but it is lifted up by quantum effects from fermionic 
zero modes [S]. Thus the model is well defined without any cutoff [2"H|26| . 

In the Lorentzian model, the bosonic action can be decomposed into two terms 

S h cx ti{F au F^) = -2tr( J F 0l ) 2 +tr(F y -) 2 (14) 

with opposite signs. The model looks extremely unstable, and hence no one has ever 
dared to study this model seriously. 

Another important difference between the Euclidean and Lorentzian models is 
the Pfaffian PiA4(A) obtained by integrating out fermions. In the Euclidean model, 
the Pfaffian is complex, and the phase plays a crucial role in SSB of SO(10) [T3UT5] . 
But it also makes Monte Carlo studies extremely difficult [I6HI8] due to the so- 
called sign problem. In the Lorentzian model, the Pfaffian is real. This is a good 
news for Monte Carlo studies, but we also lose a source of SSB. 

In the case of Euclidean model, the partition function was defined by 

Z E = J dAd^e^ s = J dAe~ Sb PiM(A) . (15) 

In the case of Lorentzian model, we define the partition function of the Lorentzian 
model by [1] 

Z L = J dAd^e tS = J dAe tSb PfM(A) . (16) 

Here we have replaced the "Boltzmann weight" e~ s used in the Euclidean model by 
e lS . This is theoretically motivated from the connection to the worldsheet theory 
When we make an inverse Wick rotation, we need to change the worldsheet 
coordinate £o = — 1^2 as well as the target-space coordinates. Applying the mapping 
rule from fields to matrices [2], we obtain (fT6|) . 

Note that e lSb in the partition function ([TB"!) is a phase factor just as in the path- 
integral formulation of quantum field theories in Minkowski space. This may give rise 
to the sign problem when one tries to study the model by Monte Carlo simulation. 
In the present case, however, the sign problem can actually be circumvented in 
the following way. The crucial point here is that the action of the type IIB matrix 
model is homogeneous with respect to the matrices. Under the scale transformation 
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i— > pA^, each part of the partition function (|16p transforms as 

5 b -> P A S h , (17) 
dA^p w ^ 2 -^dA , (18) 
PfM(A) i ^ p^^PfA^A) . (19) 

Integrating out the scale factor of the bosonic matrices first one essentially converts 
the phase factor e b into a constraint St, ~ 0. (Such a constraint is analogous to the 
one that appeared in the model inspired by space-time uncertainty principle [27..) 

It turns out that the integration over A^ in (I16[) is divergent, and we need to 
introduce two constraints 

^tr(^o) 2 <^tr(^) 2 , (20) 
l t r(A 4 ) 2 <L 2 . (21) 

This is in striking contrast to the Euclidean model, in which the partition function 
is shown to be finite without any regularization [24U26] . 

Without loss of generality, we set L = 1 in ([2T]) , and thus we arrive at the model 



Z = J ^dA^^tr(F^n)Pi^(^)^(^tr(^) 2 -l^^-ltr(A )^ , 

(22) 

where 6(x) is the Heaviside step function. Since the Pfaffian PiA4(A) is real in the 
present Lorentzian case, the model (|2"21 can be studied by Monte Carlo simulation 
without the sign problem^ Note that this is usually not the case for quantum field 
theories in Minkowski space. 



4. Monte Carlo results for the Lorentzian matrix model 

We perform Monte Carlo simulation^ of the model (1221) by using the RHMC algo- 
rithm [29] . 

It turned out that not only space but also time emerges dynamically in the 
Lorentzian matrix model. We found that the eigenvalue distribution of Ao extends 
in the large- N limit. Here, supersymmetry of the model plays a crucial role. If we 
omit fermions, the eigenvalue distribution has a finite extent, and the cutoff (|20[) in 
the temporal direction is actually not needed |30j . 

In order to extract the "time evolution" , we diagonalize Aq , and define the 
eigenvectors \t a ) corresponding to the eigenvalues t a of Ao (a = ,N) with 

the specific order t\ < • •• < ijv- The spatial matrix in this basis (i a |^i|ib) is not 



c The same procedure was also used in Ref. 1241 for simulating the Euclidean model. 
d Strictly speaking, the Pfaffian can flip its sign, but we find that the configurations with positive 
Pfaffian dominates as N is increased. Hence, we just take the absolute value of the Pfaffian in 
actual simulation. 

e For a recent review on various applications of Monte Carlo simulation to string theory, see Ref. 1281 
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Fig. 2. (Left) The extent of space R(t) 2 with N = 16 and n = 4 is plotted as a function of t for 
five values of re. The peak at t = starts to grow at some critical re. (Right) The 9 eigenvalues of 
Tij(t) with N = 16 and n = 4 are plotted as a function of t for re = 4.0. After the critical time t c , 
3 eigenvalues become larger, suggesting that the SO (9) symmetry is spontaneously broken down 
to SO(3). 

diagonal, but it turns out that the off-diagonal elements decrease rapidly as one 
goes away from a diagonal element. This motivates us to define n x n matrices 

^ b \t) = {t v+a \Ai\t v+b ) (23) 

with 1 < a,b < n and t = — J2a=i tv+a for v = 0, • • • , (JV — n). These matrices 
represent the 9d space structure at fixed time t. (This point of view can be justified 
in the large- N limit, in which more and more eigenvalues of Aq appear around some 
value t within a fixed interval St.) The block size n should be large enough to include 
non-negligible off-diagonal elements. In Fig. [5] (Left) we plot the extent of space 

R(t) 2 = -trMt) 2 (24) 
n 

for N = 16 and n = 4. Since the result is symmetric under the time reflection 
t —> — t as a consequence of the symmetry Aq —> ~Aq, we only show the results for 
t < 0. There is a critical n, beyond which the peak at t = starts to grow. 

Next we study the SSB of the SO(9) symmetry. As an order parameter, we define 
the 9x9 (positive semi-definite) real symmetric tensor 

T y (t) = itr{^ i (t)Aj(t)} . (25) 

The 9 eigenvalues of are plotted against t in Fig. [2] (Right) for n = 4.0. We 

find that 3 largest eigenvalues of T,j(i) start to grow at the critical time t c , which 
suggests that the SO(9) symmetry is spontaneously broken down to SO(3) after t c . 
Note that R{t) 2 is given by the sum of 9 eigenvalues of Tjj(t). 

It turned out that one can remove the infrared cutoffs k and L in the large- N 
limit in such a way that R(t) scales. This can be done in two steps, (i) First we 
send k to co with N as n = (3 N p (p ~ |) [3D]. The scaling behavior is clearly seen 
in Fig. [3] (Left). The scaling curve of R(t) one obtains in this way depends on (3. 
(ii) Next we send (3 to oo with L. The two limits correspond to the continuum limit 
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4 



(t-y/R(g (t-t c )/R(t c ) 



Fig. 3. (Left) The extent of space R(t) 2 for re = /37V 1 / 4 is plotted for N = 8, 12, 16 with /3 = 2. 
We plot the results against the shifted time t — t c in units of the size of the universe R(t c ) at the 
critical time. (Right) Similar plot for fixed N = 16 with re = 2, 4, 8. 

and the infinite volume limit, respectively, in quantum field theory. Thus the two 
constraints (|2U| . (|2ip can be removed in the large- N limit, and the resulting theory 
has no parameter other than one scale parameter. 

Let us discuss the second limit (ii) in more detail. We find that the inequality 
(|2"Tj) is actually saturated for the dominant configurations. Therefore, one only has 
to make the rescaling i— > LA^ in order to translate the configurations in the 
model (|2"2"|) as those in the original partition function. It turns out that R(t) for 
the rescaled configurations scales in /3 by tuning L and shifting t appropriately. In 
order to see this, it is convenient to choose L so that R(t) at the critical time t = t c 
becomes unity, and to shift t so that the critical time comes to the origin. Then R(t) 
with increasing /3 extends in t in such a way that the results at smaller \t\ scale. 
This is demonstrated in Fig. [3] (Right), where we find a reasonable scaling behavior 
for N = 16 with k = 2.0,4.0,8.0. Note, in particular, that the extent of "time" 
increases as k is increased, which is not the case in the bosonic model |30) . As we 
mentioned earlier, supersymmetry plays a crucial role for the "emergent time" in 
the Lorentzian matrix model. 

The SSB of SO(9) looks mysterious at first sight, but we can actually understand 
the mechanism quite intuitively. Let us consider the case in which k is large. Then 
the first term of (|14|) becomes a large negative value, and therefore the second term 
has to become large in order to make (fl4|) zero as required in ([22]) . Due to the 
constraint -^-tr (Ai) 2 = 1, however, it is more efficient to maximize the second term 
of (|14l) at some fixed time. The system actually chooses the middle point t = 0, 
where the suppression on Ai from the first term of (|14[) becomes the least. This 
explains why the peak of R(t) at t = grows as we increase n. 

Let us then consider a simplified question: what is the configuration of Ai which 
gives the maximum ■^■tr(Fy ) 2 with fixed tr (Ai) 2 = 1. Using the Lagrange mul- 
tiplier A, we maximize the function G = tr (Fij) 2 — Atr (Ai) 2 . Taking the derivative 
with respect to Ai, we obtain 2 [Aj, [Aj,Ai]] — XAi = 0. This equation can be solved 
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if Ai = xLi for i < d, and Ai = for d < i < 9, where are the representa- 
tion matrices of a compact semi-simple Lie algebra with d generators. Clearly d 
should be less than or equal to 9. It turns out that the maximum of -^tr (i^j) 2 
is achieved for the SU(2) algebra, which has d = 3, with Li being the direct sum 
of the spin-i representation and (N — 2) copies of the trivial representation. This 
implies the SSB of SO (9) down to SO (3). The SSB can thus be understood as a 
classical effect in the k — > oo limit. When we tune k with increasing TV as described 
above, quantum effects become important. We have confirmed |30j that the n x n 
matrix Q = ^2 i=1 A-i(t) 2 has quite a continuous eigenvalue distribution, which im- 
plies that the space is not like a two-dimensional sphere as one might suspect from 
the classical picture. 

5. Summary 

We have studied the type IIB matrix model as a nonperturbative formulation 
of superstring theory in (9+l)-dimensions. Unlike previous works, we made the 
Lorentzian matrix model well-defined by introducing infrared cutoffs instead of 
making a Wick rotation. 

Monte Carlo studies of the Lorentzian matrix model revealed the following non- 
trivial facts: 

• The infrared cutoffs can be removed in the large- N limit. 

• The theory thus obtained has no dimensionless parameter, supporting the validity 
of the model as a nonperturbative formulation of superstring theory. 

• "Time" emerges dynamically thanks to supersymmetry. 

• "Time evolution" of 9d space emerges. 

• 3 out of 9 spatial directions start to expand at some "critical time" , after which 
the SO (9) symmetry of the space is broken spontaneously down to SO (3). 

All these results for the Lorentzian matrix model suggest that (3+l)-dimensional 
expanding universe emerges dynamically from superstring theory if the theory is 
treated nonperturbatively. This may be contrasted with the quantum cosmology 
in the early 80s |32rt34] that aimed at describing the birth of the universe within 
the mini-superspace approximation^ Note also that the picture suggested here is 
quite different from that in (perturbative) superstring theory, where space-time 
with various dimensions can be obtained by compactification or by using D-brane 
backgrounds. 

The rapid expansion of the three-dimensional space observed in Monte Carlo 
simulation may be interpreted as the beginning of inflation. It would be interesting 
to investigate the microscopic origin of the inflation along this line. 

f More recently, a nonperturbative approach to quantum gravity has been pursued using the causal 
dynamical triangulation [311 - For earlier works that put forward the idea to use matrices for 
cosmology, sec Rcfs. 35, 36 See also Refs. 37, 38 for related works on emergent gravity. 
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The mechanism of the SSB relies crucially on noncommutativity of the space- 
time represented by 10 bosonic matrices. Therefore, an important issue is whether 
the usual commutative space-time appears at later times. We addressed this issue 
by studying the classical equations of motion, which are expected to be valid at late 
times [391140] . There are actually infinitely many solutions representing commuta- 
tive (3+l)-dimensional space-time. Moreover, we found a simple solution with an 
expanding behavior, which naturally solves the cosmological constant problem [ID] . 
We consider that there exists a unique solution of this kind that dominates the 
partition function of the matrix model at late times. By pursuing this direction 
further, it would be possible to understand the origin of dark energy found in the 
present cosmological observations and to predict the fate of our Universe. 

Since superstring theory is not only a theory of quantum gravity but also a 
theory of all the matters and the fundamental interactions among them, it would be 
interesting to see how the Standard Model appears at late times in the Lorentzian 
matrix model. Finding solutions to the classical equation of motions 41-45] and 
performing perturbative expansion around them would be an important direction 
as an approach complementary to Monte Carlo simulation. 

Recently, we have proposed an explicit procedure to identify the local fields corre- 
sponding to the massless modes that appear at late times [46] . The basic assumption 
is that the low-lying spectrum is essentially determined by the Nambu-Goldstone 
modes (and their extension) associated with the SSB of the (9+l)-dimensional 
Poincare symmetry and supersymmetry. The local field theory obtained in this 
way below the Planck scale has interesting generic features. It is a grand unified 
theory with the gauge group SU(fe), and all the matter fields are in the adjoint 
representation. The grand unified theory with k = 8, for instance, can accom- 
modate all the Standard Model particles. As the space-time dimensionality seems 
to be uniquely determined by nonperturbative dynamics of superstring theory, we 
consider it quite conceivable that the Standard Model emerges uniquely from this 
top-down approach. Further investigations are clearly worth while. 

To conclude, we believe that the Lorentzian matrix model, as a correct and 
tractable nonperturbative formulation of superstring theory, provides totally new 
perspectives in both particle physics and cosmology. 
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